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ABSTRACT 


In  this  work  we  investigate  the  effect  of  using  different  values  of  the  radius  in  a  constant 
modulus  algorithm  (CMA)  equalizer  for  shaped  offset  quadrature-phase  shift  keying  (SOQPSK), 
which  is  transmitted  using  the  iNET  data  packet  structure.  The  iNET  packet  structure  contains 
known  data  bits  (the  preamble  and  asynchronous  marker  (ASM)  bits)  within  each  data  packet, 
and  these  are  used  to  initialize  a  block  processing  CMA  algorithm.  Two  real-time  processing 
methods  are  investigated  here,  which  are  referred  to  as  the  ‘Single  block  processing’  and  the 
‘Block-by-block  processing’  methods.  Both  processing  methods  require  that  a  number  of  data 
packets  be  discarded  while  the  CMA  equalizer  converges.  In  this  work  we  will  investigate  the 
effect  on  the  performance  of  the  CMA  equalizer  of  changing  its  equalizing  radius,  and  also  the 
effect  of  processing  different  numbers  of  data  packets  based  on  the  same  initialization.  Based  on 
the  results  of  these  tests,  a  final  choice  of  parameters  will  be  made  and  we  will  produce  a  bit 
error  rate  curve  for  the  chosen  telemetry  channel. 

1.  INTRODUCTION 

We  investigate  the  performance  of  the  constant  modulus  algorithm  (CMA)  equalizer,  when  using 
different  values  of  the  CMA  radius,  for  which  has  data  transmitted  using  shaped  offset 
quadrature-phase  shift  keying  (SOQPSK)  modulation.  Two  different  block  processing  methods 
are  considered  for  real-time  CMA  equalization,  and  these  methods  are  applied  to  equalize  data 
bits,  which  have  been  transmitted  over  an  aeronautical  channel,  using  SOQPSK  modulation  in 
the  INET  packet  structure.  The  iNET  packet  structure  contains  known  data  bits  (the  preamble 
and  asynchronous  marker  (ASM)  bits)  within  each  data  packet,  and  these  are  used  to  initialize 
the  block  processing  CMA  algorithms.  CMA  is  a  blind  equalizer  which  recovers  unknown 
transmitted  data  by  forcing  the  equalized  points  to  lie  on  a  circle  of  a  fixed  radius  (see  [1]).  In 
this  work  we  consider  the  effect  of  varying  the  CMA  radius  on  the  performance  of  this  equalizer. 
SOQPSK  modulation  is  a  type  of  continuous  phase  modulation  (CPM)  whose  symbols  are 
therefore  of  a  constant  modulus,  and  these  actually  lie  on  a  circle  of  radius  1.  Previous  work  by 
Law  [2]  had  noted  some  difficulty  in  using  the  CMA  equalizer  for  this  application,  but  in  this 
work  we  evaluate  two  block  processing  methods,  which  are  referred  to  as  block-by-block  and 
single  block  processing  [3].  Both  of  these  methods  have  an  initialization  stage  which  requires  a 
number  of  data  packets  to  be  discarded  during  convergence  of  CMA.  We  will  investigate  the 
effect  of  using  different  values  of  radius  on  the  convergence  of  the  CMA  cost  function.  The  two 


block  processing  methods  discussed  here  have  been  chosen  in  order  to  achieve  real  time 
performanee  and  both  are  evaluated  using  the  different  radii. 

2.  COMMUNICATION  MODEL 

The  eommunication  model  for  our  applieation  is  shown  in  figure  1.  It  shows  the  transmitter, 
receiver  and  channel,  together  with  some  basic  operations  which  occur  at  the  transmitter  and  at 
the  reeeiver.  A  binary  data  stream  is  transmitted  in  the  INET  packet  structure  whieh  is  made  up 
preamble  (128  predefined  bits),  ASM  (64  bits)  and  the  actual  data  bits  composed  of  6144 
randomly  generated  binary  bits.  The  preamble  and  ASM  is  a  training  sequenee,  which  can  be 
used  to  identify  the  start  of  the  data,  to  calculate  the  MMSE  initializer  for  CMA  equalization,  and 
for  error  eorrection  [4]. 
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Figure  1 :  SOQPSK  Communication  System 

As  shown  in  Figure  1,  the  data  is  encoded  using  a  precoder,  which  converts  the  binary  input  data 
to  a  ternary  output  stream  consisting  of  {-1,  0,  +1 }.  The  signal  is  then  modulated  using  the  CPM 
modulator  with  a  pulse  shape  filter  applied  and  this  produces  the  SOQPSK  signal  which  is 
transmitted  over  the  ehannel.  At  the  reeeiver  inter-symbol  interferenee  (ISI)  due  to  the  ehannel 
effects  is  removed  by  equalization  (CMA).  Two  block  processing  methods  are  discussed  for 
processing  the  reeeived  data  while  updating  the  CMA  equalizer.  The  equalizer  output  is  then 
convolved  with  a  detection  filter,  which  acts  like  a  matched  filter  and  converts  the  equalized 
SOQPSK  signal  into  QPSK  form. 


3.  CMA  EQUALIZATION 

At  the  receiver,  the  CMA  cost  function  is  used  to  specify  a  blind  equalizer  which  is  based 
on  the  transmitted  source  signal.  The  cost  function  for  our  equalization  seheme  is  the  CMA 
eost  funetion  given  in  (1).  The  CMA  eost  funetion,/cMy4(y)’  minimizes  the  distance  of  the 
equalizer  output  from  a  circle  of  known  radius,  which  is  normally  determined  by  the  known 
eonstellation  points.  The  CMA  cost  function  is  given  by 

JcMA(.y(.n))^El^y{nf-R^y  (1) 

where  R  denotes  the  square  value  of  the  CMA  equalizer  radius  and  y(n)  represents  the  equalizer 
output  block  at  time  n.  It  is  known  that  SOQPSK  modulation  symbols  of  the  transmitted 
modulated  data  are  of  a  constant  modulus  and  of  radius  1,  but  in  this  work  four  different  values 
of  this  CMA  radius  squared,  (R=l,2,3  and  4)  are  evaluated  for  equalization  performance. 

For  CMA  radius  i?=l,  figure  2(a)  and  2(b)  shows  respectively  the  equalized  output  and  the  QPSK 


constellation  after  the  detection  filter  is  applied  at  the  receiver  at  SNR  =18  dB  for  channel  1. 
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Figure  2(a):  Equalized  SOQPSK  for  radius  1 

When  a  CMA  radius  squared  value  of  /?  =2  i 
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(b)  Equalized  QPSK  constellation  for  radius  1 
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Figure  3:  (a)  Equalized  output  for  radius  2  (b)  Equalized  QPSK  constellation  for  radius  2 

The  symbols  of  the  equalized  SOQPSK  modulated  symbols  appear  to  lie  on  a  circle  of  radius  a/2 
while  the  QPSK  constellation  points  are  located  at  +!-  V2.  Similar  results  have  been  observed  for 
the  higher  values  of  radius,  but  for  the  radius  squared  values  of  /?  =  3  and  4,  the  SOQPSK  circle 
radius  and  constellation  points  actually  lie  at  the  values  of  3  and  4,  respectively. 

4.  IMPLEMENTATION  OE  CMA  EOR  AERONAUTICAL  TELEMETRY 

In  order  to  evaluate  the  equalizer  performance,  we  will  use  the  previously  measured  aeronautical 
telemetry  channel  1  for  simulation  and  testing  [5].  This  channel  is  of  length  9  and  was  measured 
on  the  taxiway  before  the  aircraft  takes  off.  Flowcharts  are  presented  below  for  the  two  real-time 
block  processing  methods  implemented  in  this  work.  Section  4.1  describes  block-by-block 
processing  and  section  4.2  describes  single  block  processing. 


4.1  Block-by-Block  (BxB)  Processing  Method 


Figure  4:  Flow  chart  of  BxB  processing 

For  BxB  processing  the  CMA  equalizer  is  updated  iteratively  on  multiple  packets.  Errors  are 
counted  after  a  specified  number  of  packets  have  been  processed  to  allow  for  convergence,  and 
averaging  is  done  over  a  total  number  of  specified  packets.  The  flow  chart  in  figure  4  shows  the 
implementation  of  the  block-by-block  processing  method. 


4.2  Single  Block  Processing  Method 
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Figure  5:  Flow  chart  of  Single  block  processing 

For  single  block  processing  the  first  packet  is  used  to  update  the  equalizer  filter  coefficients  until 
convergence  of  the  equalizer.  The  error  analysis  is  then  done  by  bringing  in  new  packets  to 
update  the  equalizer. 


5.  PERFORMANCE  EVALUATION  RESULTS 


In  this  section  we  present  performance  evaluation  results  for  the  CMA  equalizer  for  channel  1. 
First  the  convergence  curves  for  the  block-by-block  processing  and  single  block  processing 
methods  are  compared  for  different  values  of  CMA  radius  in  order  to  determine  the  number  of 
packets  required  for  initialization  of  each  block  processing  method,  and  then  we  compare  the  bit 
error  rate  performance  when  using  different  values  of  CMA  radius. 

5.1  Evaluation  of  different  values  of  CMA  radius 

In  this  section  the  values  of  CMA  radius  is  changed  and  we  compare  the  performance  of  the 
CMA  equalizer  by  plotting  the  average  bit  error  rate.  First  we  consider  the  CMA  cost  function, 
Jcma  to  determine  the  number  of  data  packets  required  for  convergence  of  the  equalizer.  Figures  6 
(a)  and  (b)  show  for  a  single  Monte  Carlo  run  the  convergence  curves  of  Jcma  for  BxB  processing 
and  for  Single  block  processing  respectively  for  radius  squared  values  of  R=l,  2,  3  and  4. 
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Figure  6:  Convergence  curves  for  Jcma  using  radius,  i?=l,  2,  3,  and  4 

The  number  of  data  packets,  which  must  be  discarded  during  initialization  for  these  two  real¬ 
time  methods,  can  be  determined  from  the  convergence  curves  of  figures  6  (a)  and  (b).  We  note 
that  the  cost  function  value  increases  with  increasing  radius.  In  addition  the  larger  the  value  of 
the  radius  the  more  iterations  it  takes  for  both  methods  to  converge.  50  packets  as  startup  (for 
convergence)  appears  to  be  sufficient  when  using  a  radius  squared  value  of  1  and  2  for  this 
particular  Monte  Carlo  run,  while  radius  3  requires  about  150  packets  for  convergence,  and 
radius  4  requires  about  250  packets  for  convergence  for  both  methods.  The  number  of  packet 
iterations  needed  for  initialization  of  the  block-by-block  processing  method  is  comparable  to  the 
number  of  iterations  required  for  the  single  block  method.  Table  1  summarizes  the  parameters 
that  were  chosen  from  the  convergence  curves  of  figure  6,  in  order  to  perform  the  bit  error 
analysis  for  the  four  different  values  of  CMA  radius.  These  results  presented  in  figure  7  use  the 


last  column  of  Table  1  so  that  the  number  of  packets  used  for  error  averaging  is  consistent  for  the 
two  methods. 


Table  1:  Number  of  Packets  used  for  Radius  Performance  Analysis 


Radius 

Total  Block  Length 
(BxB) 

Convergence 

iterations 

Average  errors  over 

1 

700 

50 

650 

2 

700 

50 

650 

3 

750 

150 

600 

4 

850 

250 

600 

In  figure  7  the  bit  error  rate  curves  for  the  two  methods  are  plotted  for  channel  1  based  on  the 
packet  initialization  specified  in  Table  1  and  determined  from  figures  6.  The  theoretical  lower 
bound  curve  for  the  AWGN  channel  is  included  to  show  the  best  achievable  equalizer 
performance  for  a  single  tap  channel  with  additive  white  Gaussian  noise.  Note  that  the  best 
results  are  produced  for  CMA  radius,  R  =  2.  A  baseline  curve  for  the  performance  of  the  MMSE 
CMA  equalizer  is  also  included.  This  curve  was  produced  by  taking  a  Monte  Carlo  average  of 
the  final  error  rate  of  an  MMSE-initialized  CMA  equalizer  where  convergence  is  achieved  using 
a  single  fixed  packet  [4].  It  is  not  a  real-time  procedure  and  it  produces  almost  identical  results 
for  both  radius  squared,  R  =  I  and  2. 
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Figure  7:  BER  vs.  SNR  for  different  radius  values  (R=  1,  2,  3,  and  4)  for  channel  1. 

We  note  from  figure  7  that  for  both  methods  the  radius  i?  =  1  curves  produce  the  highest  bit  error 
rates  at  low  SNR  followed  by  radius  squared,  R  =  2,  and  both  achieve  the  baseline  performance 
of  the  MMSE  CMA  equalizer  described  above.  At  low  SNRs,  the  radius,  R  =  2  and  4  curves 
produce  lower  BER  than  even  the  baseline  MMSE  CMA,  but  at  higher  SNR  their  error  rates  are 
less  stable  and  are  higher  than  the  baseline  MMSE  CMA.  In  addition  radius  /?  =  3  and  4  require 


higher  numbers  of  packets  to  be  discarded  during  initialization  when  compared  to  radius,  i?  =  1 
and  2.  For  the  reasons  stated  above,  we  choose  to  continue  the  performance  evaluation  of  CMA 
by  comparing  only  radius,  R  =  I  and  2. 

5.2  Block  Length  Evaluation 

In  this  experiment,  we  test  the  effect  of  increasing  the  block  length  in  BxB  processing  method 
for  radius  1  and  radius  2  using  channel  1,  by  plotting  bit  error  rate  curves  for  different  packet 
lengths.  These  results  are  shown  in  figure  8. 
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Figure  8:  Block  length  comparison  for  Block-by-Block:  Radius  1  and  Radius  2 

CMA  squared  radius,  R  =  2  curves  provide  lower  bit  error  rates  than  the  radius  1  curve  and  the 
performance  using  R  =  2  remains  consistently  stable  as  the  packet  length  varies  from  500  to  2000 
with  improved  performance  at  higher  SNRs  for  a  packet  length  of  2000.  The  radius  1  curves  are 
more  varied  for  different  packet  lengths,  converging  to  a  smoother  curve  for  a  packet  length  of 
2000.  The  Single  block  processing  method  shows  performance  similar  to  this. 

5.3  Comparing  BxB  with  Single  Block  processing  methods 

In  this  experiment  we  compare  the  block-by-block  processing  method  with  the  single  block 
processing  method  for  the  aeronautical  telemetry  channel  1.  Both  block  processing  methods 
(with  radius,  i?  =  1  and  2)  will  use  an  initial  set  of  50  packets  for  startup  (and  convergence)  of  the 
equalizer,  so  these  packets  are  discarded.  Errors  are  then  counted  and  averaged  over  a  total  of 
650  data  packets  for  a  total  of  700  data  packets.  Figure  9  shows  the  error  rates  for  the  two  real¬ 
time  methods  discussed  here. 


Channel  1 


SNR  (dB) 

Figure  9:  BER  vs  SNR  for  Channel  1 

We  observe  that  CMA  radius  squared,  R  =  2  produces  consistently  better  results  for  this  channel. 
The  performance  of  BxB  method  is  very  similar  to  that  of  single  block  method  for  radius,  R=2, 
but  there  is  more  of  a  difference  in  the  performance  of  the  two  methods  when  CMA  radius  1  is 
used. 


6.  CONCLUSIONS 

In  this  work,  we  investigated  the  effect  of  using  different  values  of  radius  squared,  R  for  a  CMA 
block  equalizer  for  SOQPSK,  which  has  been  transmitted  using  the  iNET  data  packet  structure. 
The  known  data  in  the  iNET  packet  structure  is  used  to  determine  the  MMSE  equalizer  which  is 
used  to  initialize  the  CMA  equalizer.  Two  implementations  of  the  block  processing  CMA 
equalizer  have  been  investigated,  and  these  have  been  shown  to  have  similar  performance.  Both 
of  these  methods  are  implemented  to  achieve  real-time  equalization  of  the  transmitted  data,  and 
they  both  require  that  a  number  of  data  packets  are  discarded  during  convergence  of  the  CMA 
equalizer. 

In  general,  a  CMA  radius  squared  value  of  R  =  2  would  be  our  best  choice  for  the  CMA 
radius,  since  it  achieves  consistently  lower  bit  error  rates  over  the  full  range  of  SNRs  than  radius 
1.  Eor  the  radius  squared  values  of  R  =  3  and  4,  the  bit  error  rates  at  low  SNR  can  be  lower  than 
that  of  radius  squared  R  =2,  but  these  require  a  larger  number  of  data  packets  to  be  discarded  for 
initialization  and  the  high  SNR  behavior  for  these  higher  radius  values  appears  to  be  unstable 
and  it  diverges  from  the  baseline  performance  provided  by  the  MMSE  CMA  equalizer  curve. 
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